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Abstract 



It is generally agreed that interstellar dust grains consist of two main components, 
namely, silicates and graphites. Some models, like MRN model, assume these grains to 
be homogeneous spheres following a power law size distribution. This paper presents, 
in the framework of Mic theory, a parametrization of extinction spectrum curves of 
the silicates and the graphites separately in terms of frequency and the minimum 
and maximum of sizes in the distribution. Analytic expressions in ultraviolet and far- 
ultraviolet are presented for both types of grains. The values of maximum and minimum 
sizes for which these equations are valid have been identified. These equations can be 
useful in a number of situations involving silicate and graphite grains. 
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1 Introduction 



In a recent publication we analyzed extinction spectrum of a collection of homogeneous 
spherical particles of a unknown size distribution [1] ( hereafter referred as RS-I). It was 
shown that an extinction spectrum, in general, has some easily identifiable characteristic re- 
gions where the extinction-frequency relationship can be approximated by simple empirical 
formulae involving first four moments of the particle size distribution. The analysis clearly 
exhibited the manner in which essential features of the particle size distribution gets coded 
into its extinction spectrum (more generally known as the Interstellar Extinction Curve in 
the astrophysical situation). It was demonstrated that the moments could indeed be ob- 
tained from its extinction spectrum and that with the knowledge of these moments it was 
possible to reconstruct the size distribution function. 

The analysis in RS-I however, assumed that all particles were of same material for which 
the refractive index did not vary with frequency. The question arises what happens when 
a collection of particles has more than one material component with their refractive indices 
varying with wavelength. We address ourselves to these questions in this paper. In other 
words, the purpose of this work is to examine the possibility of extending the ideas developed 
in RS-I to a multi-component system where refractive index of the individual components 
varies with wavelength of incident electromagnetic radiation. 

One suitable model where these extensions of RS-I can be studied is a modified version 
of interstellar dust model due to Mathis, Rumpl and Nordsieck (MRN) [2]. The MRN model 
assumes that interstellar dust consists of individual grains of homogeneous spheres of sili- 
cates and graphites having very definite power-law size distributions. For graphite grains, 
there is a "| — |" approximation [3]. This modification essentially converts two compo- 
nent MRN model to a three component MRN model. Investigations in this work show that 
the frequency-size relationship can still be obtained even when the refractive index of the 
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grains is varying with wavelength. But, because of rapid variation of refractive index with 
frequency, now one needs to divide the spectrum in suitable intervals and obtain separate 
relations for each sub-interval. This paper presents these relationships in ultraviolet and 
far-ultraviolet frequency domains for silicate and graphite grains individually. Thus, we not 
only demonstrate, the successful extension of ideas in RS-I to materials whose refractive 
index varies with wavelength, but in the process we also present relationships which can be 
used for testing various models of interstellar dust. 

This paper is organized as follows. Section 2 describes the dust model considered here. 
Section 3 gives the functional form of extinction in terms of size distribution parameters as 
well as the frequency for both graphite and silicate grains. Subsequently, numerical results of 
extinction obtained from these formulae are compared with exact Mie theory computations 
in section 3. The possible use of such formulae for individual components to multi-component 
system has been demonstrated in section 4. Finally, we conclude by summarizing and dis- 
cussing the results of this paper in section 5. 



A classic model of interstellar dust was proposed more than 30 years ago by Mathis, Rumpl 
and Nordsieck (MRN) (1977). Since then, the basic model is being used even today albeit 
with some modifications. It has not been fully superseded by later studies. The MRN 
model uses two separate populations of bare silicate and graphite grains with a power-law 
distribution of sizes of the form: 



where a denotes the radius of the spherical grain varying within the limits Gq (minimum 
radius) and (maximum radius). The size range for graphite and silicate grains is: 



2 The Dust Model 
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Graphite grains : Oq ~ 0.005/xm a. 



"m 



0.25/xm, 



(2) 
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Silicate grains : oq ~ 0.005//m a. 



"tn 



0.25//m. 



(3) 



Further, the graphite material is supposed to be present in two distinct structural varieties 
within the specified range of Oq and a^- This plausibility lies in the fact that graphite 
is highly anisotropic material. The refractive index of graphite, therefore, depends on the 
orientation of electric field relative to the structural symmetry. Owing to difficulties in 
calculations of exact scattering quantities due to anisotropy, workers have taken resort to an 
approximation known as "| — |" approximation [3]. In this approximation graphite grains 
are represented as mixture of isotropic spheres, | of which have refractive index m — m\\ 
(referred as graphite parallel) and | have the refractive index m = m±_ (referred as graphite 
perpendicular) . This modification effectively makes MRN model a three component model. 

3 Extinction as a function of frequency and size distri- 
bution parameters 

Exact extinction coefficient, K^xti was obtained using the formula. 



where Qext{x) is the extinction efficiency of an individual scatterer of size parameter x and 
X = 27ra/A with A as the wavelength of the radiation. The exact extinction efficiency Qg^t 
for a spherical homogeneous scatterer can be computed using Mie formulas. The necessary 
refractive index particulars for various components (at various wavelengths) were taken from 
the tables provided by Draine [4] on his website. The number of particles per unit volume, N 
has been arbitrarily fixed at = 4.4 x 10^. in these calculations. This being a multiphcative 
constant, its admitted value does not make any effective difference in the functional form of 
Kext we wish to determine. 




(4) 



In case of size distributions of spherical particles having constant refractive index (wave- 
length independent) it was shown in RSI that in the K^xt — ^ graph, in general, we can 
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identify frequency intervals where Kext has distinctfy hnear or parabohc or asymptotic be- 
haviour having corresponding functional forms of Kext'- 

K^xt ^ Pi{m)a'^i'^ + p2{m)a^i' + p3{m)a'^, 

etc., where, li,l2,Pi,P2,P3 are coefficients which have arbitrary dependence on refractive in- 
dex m. a^, and a"' are the 2nd, 3rd and 4th raw moments of the size distribution /(a). 



In the present situation, where m varies with the coefficients Zi, /2,Pi,P2,P3 will also 
show variation with u and hence the distinctly simple linear and parabolic forms of K^xt will 
not be observed within a meaningful frequency interval. Hence, our approach has been to 
study the K^xt — ^ graph for each material component separately for various ranges of a, 
which means varying oq, a^. In case of the power-law size distribution, /(a) = ca~^'^, we 
have relations: 

5 5/2/, 1 \-l 

5/ 1\/ 
a2 = 5a^(l-i)(l- 



1 \-l 



^ TV' 



= hal{n-l){l ^ ^ ^ 



where n — {am/ a-oY^'^ ■ Consequently, with a functional form Kext{oio-,o.m-,^)-, our investiga- 
tions reveal that the extinction in the UV and FUV regions for the materials considered have 
the following general form: 



ao, v) + il){am, y) ■ (5) 



The functions (f) and i}) have forms which change in various frequency sub-intervals. For 
each component of graphite (parallel and perpendicular), four formulae were needed to fit 



extinction. For silicate, three formulae were sufficient. Several values of and oq were 
considered. Formulae presented here are valid in the size limits: 



Observations of the spectra suggested that the functions 0(ao, ly) and ^(a^, v) can have the 
simple forms: 



where oq, Om are taken in units of 10~^cm, u in units of lO^cm"^. The number of significant 
terms contributing to extinction in and ^ depend on the material as well as the frequency 
interval considered. The forms of 0(ao,i^) and ^(a^,i/) have been constructed by careful 
analysis of the regional (frequency sub intervals) behaviour of the extinction for each of the 
materials considered. In the process of developing analytic formulae, care has been taken to 
have a good compromise between accuracy and calculational simplicity so that our analysis 
could be applied expediently for purposeful dust modelling within the power law framework 
considered here. 

Following relationships have been obtained: 

3.1 Homogeneous graphite grains with refractive index m = m± 

1. For 1000 < A < 1460A (FUV region I): 



Graphite grains : 0.002/xm < Oq < 0.005/im; 0.15/xm < a^ < 0.25/xm, 



(6) 



Silicate grains : 0.004/xm < ao < 0.006//m; 0.2/xm < < 0.4/xm. 



(7) 





K = Cal^'^ (0.259 + 20.3073(i/ - 0.8485)^) - (i/ao)^/^(o.267z/ - 0.16048+ 




14.8274(1/ - 0.8428)^(1/ - 0.67905)) - 0.0488(^ + 



0.1551 



) 



(8) 
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2. For 1460 < A < 1900A (FUV region II): 



K = Ca 



5/2 



1.97754 



2.42703 0.79483 a 



1/2 



- + - (0.2164 - 0. 



278121/ + 



0.43347 



:i.o- 



0.6856 



) 



0.5263, X ao /, 0.606,2 0.03173, 0.6262, 

1.0 ) ^(1.0 f + 1 46.52t/ 1/ - 0.6856 

V ' ^ V V V 



[y - 0.5263)(i/ - 0.606)) - 0.0488(^ + ^^^^^)(0.8667 + 24.1(i/ - 0.6112)^) 

' am 



(9) 



3. For 1900 < A < 2500A (UV region I): 



5/2 



+ 100|i/- 0.4654|=^-i/ay^(- 



1 



+ 



2.7 + 490(i/ - 0.4654)2 

11307(i/-0.4673)^-258.26i/(i/-0.4759)2) -ao f . "'""^'^ ■ A 

^ ' ^ ^' Vi.o + 6.9(1.0 - ^^) + 76.187i/|i/- 0.45661) ^ 



1058 + 6.5985|z/ - 0.47361 
0.0813 



1 2398 

- 0.0551l(^ + ) (1.9503 - X.lVlv + 0.583|i/ - 0.452 

4. For 2500 < A < 4000A (UV region II): 



(10) 



z/^/2(o.ii61 + 0.60i/ - 1.9684i/^/^(0.4 - vfl'^{v - 0.25)) - 0^^^(0.461+ 



1.7567|i/ - 0.325| + 180.6(i/ - 0.325)^) - aoi^^ (20664(i/ - 0.3244)^+ 



1.30191/ 



N / 1 0.278x 

,)-0.0507(-^ + ^) 



1.0 + 304.9(1.0-^)2 
3.2 Homogeneous graphite grains with refractive index m = m 

1. For 1000 < A < 1120 A (FUV region I): 



K = Cao 



5/2 



t/2/=^ (0.3473 - 6.284(1 - i/)2) - ayV^/2(2.3047t/ - 2.0047+ 



1.38(1 - ^)(i - 1)./.) _ „„(5.55(i - 1)"=(1 - ^)|1 - 



0.2979 

1/1/2 



11.65, /0.046 O.Ollx 



(12) 



2. For 1120 < A < 1500A (FUV region II): 



K = Ca'^' 



0.06685 - 0.06056i/ + 0.2771/^ - 0^^(0.154821/ - 0.07954)- 



2 1 

ao(0.75484i/ - 0.495) - 0.018(^ + — 

am ^rn 

3. For 1500 < A < 2400A (FUV region III): 



K = Ca^^ 



1 "^/^ 

— 7;rf0.0864 + - 0.4675)^ - 7.16(i/ - 0.49)^) - ^fo.0065+ 



- 0.4922)^ - 1.76(. - 0.6112)') - a„(0.0054 + '>-n- -^-^-^<^^\ 
V ^ ' ^ ^ y "V 1.0 + 118(nM- 0.5405)2 



0.03913 + 0.07925|t/ - 0.5311| 0.008 - 0.2116(i/ - 0.5486)^ 



1/2 
dm 



var. 



4. For 2400 < A < 4000A (UV): 



K = Co; 



5/2 



0847 

0.1495 - ' ^ |t/ - 0.36751 - ao/^(i/(0.1915 - 0.817|i/ - 0.3719|)- 



33 



.0|i/ - 0.36751(0.4167 -v){v- 0.25)^) - oo — 



0.00356 



(1 + 15.017|z/- 0.4031] 



0-26 , 0.00436 + 0.0122i/|i/- 0.35421 
^^z/(1.0 - 1.36i/) ■ - 



Qml^ 



3.3 Homogeneous silicate grains 

1. For 1000 < A < 1460A (FUV region I): 



K = Cao 



5/2 



0.11668 - 0.02745t/ + 0.171/^ - 0^^(0.35421/ - 0.2110)- 



- ao(o.6677z/ - 0.3733) + (2a„(l - z/)(z/ - 0.6856) - ^ - °'°°^^ ) 
2. For 1460 < A < 2000A (FUV region II): 



K = Cao 



5/2 



0.29168i/ - 0.02268 - 0^^(0.06341/ - 0.01623)- 



/ X /0.05 0.0062 

ao(0.64634z/ - 0.3478] - [-^ + 



Ctm 



(17) 



3. For 2000 < A < 4000A (UV): 

K = CaJ^ 



-0.02833 + 0.197Z/ + ^^^^ + 0.006567a„z/ - 



0.016417/ 1-411 



0.435 I 



al 1(1.0 + 1011-^^1)2 

In all the above formulae fl8|)- f|T8|) . C = ttN/IO^; oq, flm are in units of lO^^cm; u in units of 
lO^cm"^ 

4 Numerical comparisons 

Figure 1 shows a comparison of exact extinction curves with predictions of formula for 
graphite grains. The " | — |" approximation has been used in these comparisons. The three 
graphs presented in this figure are for oq = 0.002/im, 0.0035/im and 0.005/im with fixed value 
of = 0.25/im. It can be seen that the agreement between exact results and predictions of 
formulas is excellent. Although not shown here, it has been noted that while the variation 
in oq results in bigger changes in extinction, the variation in extinction due to variation in 
Qm is comparatively much weaker. 



For silicate grains, the comparison of predictions from f|T6l) . f|T7|) and f|T8l) have been compared 
with exact results of computations from (jlj) in figure 2. In this figure, the three admitted 
values of are 0.004/im , O.OOSyum and 0.006yum with fixed am = 0.30/im. Predictions of 
formula can be seen to be extremely good. Here also the variation in extinction with am is 
much weaker in comparison to that with ao. 



To give an idea of errors in predictions of the formulas obtained, we have displayed per- 
cent error details in figures 3, 4 and 5 respectively for graphite perpendicular, graphite 
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parallel and silicate grain formulae. The percent error has been defined as: 

[Kexti^xact) — K^j-J formula)] x 100 

Percent error — — — r 

Kext{exact) 

In general, the error is within 1% but for a few regions it can go upto about 3%. The results 
are representative results for particular values of oq and a^. However, It has been varified 
that this result is valid for entire range of Oq and values considered in this work. 

5 Conclusions and discussions 

In this work, we have presented the extinction spectrum analysis for astronomical silicate 
and graphite grains in the wavelength range ~ 1000 - 4000 A which covers the FUV and 
UV regions. The grain size distribution range covered by us is in keeping with the accept- 
able compositional aspects of interstellar dust models corresponding to average interstellar 
spectra data. From this point of view, we feel that the results and the analytic formulae 
presented in this paper have application-worthiness in the sense that accurate estimates of 
extinction contributions for these material components in the FUV and UV regions can di- 
rectly be made. Needless to say that a fuller and more complete account of the extinction 
contributions made by each and every material components in respect of the various dust 
models are needed for ascertaining their effectiveness in reproducing the average interstellar 
dust extinction spectrum. This necessitates further large -scale investigative work-outs and 
mathematical analyzes of similar nature as provided by us in this work. 

However, as is mentioned in section 1, our prime motivation behind this work has been to 
analyze the extinction spectrum for a size distribution of spherical particles having variable 
refractive index (m = m(A)). Earlier, the same problem with constant m (at all wavelengths) 
was analyzed to have simple extinction-frequency relationships in different frequency ranges 
involving the first few lower order moments as well as some coefficients which seemed to be 
arbitrarily dependent on the refractive index. In case of power- law distribution /(a) = a~^^^ 
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all the moments (a, a^, a^, a^) become simple functions of the end points oq, Om- Con- 
sequently, we have looked for functional forms of extinction in terms of variables oq, Qm 
and 1/(1/ X) in various sub-regions as has been given in section 3. This analysis can also be 
extended to more complex distributions that have been used by some authors [5,6]. In the 
present work, we have restricted ourselves to UV and FUV regions and power-law distribu- 
tion for substantive demonstration of workability of the formulae developed along the line 
of our earlier work. It is our wish to extend this further to cover the visible and infrared 
regions of extinction spectrum for both these grain components as well in near future. 
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Figure captions 



Figure 1. Comparison of predictions of equations ([8])- f|TT]) and f|T2l) - f|T5l) in "| — i" approx- 
imation with exact computations from (jlj). Solid lines are predictions and points are exact 
computations. In this figure, am = 0.25/im and ao = 0.005/im (sohd line), ao = 0.0035yum 
(large dashed line) and = 0.002/xm (small dashed line). 

Figure 2. Comparison of predictions of equations f|T6l) . f|T7|) and f|T8l) with exact compu- 
tations from (jlj). Solid lines are predictions and points are exact computations. In this 
figure, am = 0.30 fim and ao = 0.006/im (solid line), ao = 0.005/im (large dashed line) and 
ao = 0.004/im (small dashed line). 

Figure 3. Percent errors in equations (IH]), (Ej), ( fTOl) and ( fTTl) with respect to exact com- 
putations from (jl]). In this figure, a^ = 0.25/im and ao = 0.005/im. 

Figure 4. Percent errors in equations (fT2|) . (fT3|) . (IT^ and (fT5|) with respect to exact compu- 
tations from (jl]). In this figure, a^ = 0.25/im and ao = 0.005/tm. 

Figure 5. Percent errors in equations f|T6|) . f|T7|) and f|T8|) with respect to exact computations 
from (jlj). In this figure, am = O.SOfim and ao = 0.005/im. 
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